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A thermowell is a metallic product fitted into the wall of a pipe or 
vessel so as to permit introduction of a thermometer or thermocouple for 
the purpose of measuring the temperature of the contents. It is designed 
so as to maintain the integrity of the pressure boundary without introducing 
unacceptable measurement errors or time lags. This monograph summarizes 
the results of analytical studies made by the writer during the past two 
years of the mechanical/structural integrity of thermowells. It is obvious 
that much of this material is also applicable to other insertions such as 
sampling tubes, which, however, need not sustain the differential pressuriza- 
tion to which thermowells are subject. 

An existing section (1) of the ASME Power Test Codes, based very 
largely upon analysis reported by J. W. Murdock (6) , represents the consensus 
of the ASME Committee PB51 (on Thermowells). Designers have recently found 
it difficult to reconcile the strict requirements of this document with the 
practical necessity of providing thermowells for boiler feed discharge and 
main steam services. In the summer of 1972, Mr. J. E. Leary, Chief Control 
and Instrumentation Engineer of Bechtel Power Corporation, asked the writer 
to examine the structural integrity of thermowells and to compose recom- 
mendations for analysis of high pressure thermowells. A report (3) and 
a supplement (4) were produced shortly thereafter. A related study (5), by 
Professor T. M. Houlihan, examining the thermal performance of thermowells 
was also produced at this time; from this study it may be inferred that 
thermowell tip details which permit full assurance of structural integrity 
impose no problems of inaccurate temperature measurement or thermal time 
lag. 

Subsequent to the production of the reports (3), (4), and (5), Mr. 

Leary asked the writer to solicit and collect comments from Mr. Murdock and 
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members of the PB51 Committee. The writer is very pleased to acknowledge 

the participation and cooperation of the following persons: 

Mr. J. W. Murdock, formerly of the Applied Physics Department, 

U. S. Naval Ship Engineering Center, Philadelphia, and presently 
a private consulting engineer. 

Mr. L. A. Dodge, Bailey Meter Company 

Mssrs. J. Archer and T. Reitz, Gilbert Associates 

Mssrs. R. F. Abrahamsen and J. D. Fishburn, Combustion Engineering, Inc. 
Mssrs. A. Lohmeier and A. J. Partington, Westinghouse Electric Corporation 
Mr. W. N. Wright, TemTex Temperature Systems and Components 
Mr. W. 0. Hays, ASME (Secretary, PB51) 

The most significant change recommended by the writer in his first 
report (3) was a drastic relaxation of the requirements with respect to 
simple pressurization. This matter is discussed in detail in Appendix A 
hereto. In the first report, in what seems to have proved to be a mistaken 
attempt at simplifying the presentation of the analysis, the writer consid- 
ered internally pressurized hollow cylinders, a case to which most pertinent 
engineering literature has addressed itself. However, the point was 
established that with any material failure theory which is independent of 
the first scalar invariant of the stress tensor, the internally and externally 
pressured cases are equivalent if strains are limited to the order of magni- 
tude of elastic strains. In the present Appendix A, the analysis is specifi- 
cally directed to the externally pressurized case. However, precisely the 
same results are obtained as previously. 

The comments generated in response to the writer's request dealt 
preponderately with the matter of pressurization and its consequences. 

While the need is acknowledged in the case of supercritical pressure instal- 
lations to depart from the limitations imposed by strict application of the 
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rules for externally pressurized vessels to be found in Section VIII of 
the ASME Boiler and Pressure Vessel Code, there is a reluctance actually 
to do so, and some responders feel that a test program is called for. 
However, the writer is absolutely certain in his own mind that internal 
pressurization tests would be absolutely useless. External pressurization 
tests might indeed prove useful as a basis for establishing rules even more 
liberal than the writer suggests (cf. Appendix A) should the need to do so 
ever arise. The writer's first report failed to cite the classic experi- 
ments of Prof. P. W. Bridgman (2) which should be adequately demonstrative 
of the ability of very thick wall, externally pressurized cylinders to 
retain pressure integrity even under extreme pressurizations. The comments 
of Mr. J. D. Fishburn conclude with the statement: 

"...In fact. Professor Brock's solution remains a conservative solu- 
tion apart from the Code limitation of S = (5/8) a. This is for 

m 

three reasons. Firstly, as stated by Professor Brock, 'for external 
pressurization. .. gross deformation is such as to modify the geometry 
advantageously,' secondly, a non-workhardening material is assumed, 
and, thirdly, there is considerable experimental evidence to show 
that the generalized Tresca condition for yield is of itself conser- 
vative. 

"For these reasons there should be no objections from the 

industry to this section of Professor Brock's analysis, particularly 

since the Nuclear Code specifies a value of 3S as the limiting 

m 

stress range (Section III NM3222.2)." 

Briefly, the writer, having studied all comments on the matter of 
pressurization and rules relating thereto, remains firm in his conviction 
that the liberalized criterion he recommends (cf. Section 10 of Appendix A) 
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is still so conservative that it might be further liberalized if the need 
to do so should ever arise. 

A second, and relatively unimportant pressure criterion relates to 
the thickness of tip closure for cantilevered thermowell. The requirement 
stated in Section II of Appendix A is simple and probably very conservative. 

The third recommended criterion relates to the necessity of assuring 
that the mechanical excitation provided by the forces exerted on the thermo- 
well by fluid flowing past it not be in resonance with a natural vibration 
mode of the thermowell structure. The analysis of this matter naturally 
divides itself into three aspects: (a) estimation of the exciting frequency, 

(b) estimation of the response (natural) frequency of thermowell vibration, 
and (c) consideration of how closely these may be permitted to approach one 
another. 

The most troublesome and controversial of these subproblems is the 

first, which is discussed in detail in Appendix B hereto. For many years 

it was thought that the dimensionless parameter (called the Strouhal number) 

determining the frequency of vortex shedding for fluid flow past a fixed, 

rigid cylinder had a definite value (N^ = .21, appr. ) for all flows faster 

than those characterized by Reynold's number (based on cylinder diameter) 

N = 800 (appr.) and Murdock's analysis (6) is based upon N = .21. How- 
R 

ever, more recent data, not available at the time of Murdock’s study, 

indicate that the Strouhal number may become as great as = .45 for 

large Reynold's numbers (N = 6.2 x 10 6 , appr.). This suggests that the 

R 

frequency of excitation can become more than twice as great as previously 
contemplated. The writer's first report (3) employed this newer information 
in estimating the excitation frequency. 

Comments received in this regard pointed out that there was lack 
of coherence in the vortex shedding which takes place for these higher 
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values of Reynold's number and that, accordingly, the danger of resonant 
excitation is thereby diminished. There can be no argument with this 
contention but the question remains: to what extent is the probability of 

resonant excitation actually reduced? The experimental data invite a varie- 
ty of interpretations and, in the writer's opinion, it is prudent to assume 
that coherent excitation of sufficient duration (perhaps as briefly as a 
fraction of a second at the frequencies involved) to cause damage can 
occur if excitation * frequency based on N = .45 equals or exceeds the 

O 

natural response frequency. Accordingly, the criteria recommended herein 
are based upon the possibility of coherent excitation with N = .45 in the 

O 

appropriate range of Reynold's numbers. However, if these criteria are 
not met, it is suggested that the designer-engineer feel free to re-examine 
the matter, making use of whatever new information may at that time be 
available and making a special examination of the probability and conse- 
quences of coherent excitation. In correspondence with the PB51 Committee, 
the writer has suggested that the use of appropriate vibration test instru- 
mentation applied to existing thermowell installations can possibly provide 
information which will help to assess the dangers of excitation in this 
regime. It should be made clear that this area is one in which persuasive 
information is indeed lacking and the writer's recommended criteria are 
intended to be definitely conservative. 

The second of these three subproblems is that of estimating the 
natural response frequencies of thermowell vibration. This matter is the 
subject of Appendix C hereof, in which, throughout, it is assumed that we 
are dealing with a thermowell which is firmly attached to the pipe or ves- 
sel wall at one end (the root) and is free at the other end (the tip) . 
Section 14 of Appendix C discusses a structural mode, not elsewhere 
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discussed except in the supplement (4) to the writer's earlier report, 
which involves ovalization of the pipe or vessel and which is at relatively 
low frequency. An argument is offered for regarding this mode as of no 
practical significance, but it would certainly be desirable to have experi- 
mental evidence that this is the case. 

The bulk of Appendix C relates to the estimation of the lowest mode 
of cantilever beam vibration. There are several complicating factors. The 
structure is nonuniform and is so short and stubby that rotatory inertia 
and shear deformation have a significant effect in reducing the response 
frequency as compared to what might be calculated by the use of so-called 

I 

"elementary" theory. Furthermore, the pipe or vessel wall to which the 
thermowell is attached is itself flexible and possesses mass and this causes! 
a further reduction in response frequency, compared to the elementary 
assumption of fixed root. 

In Appendix C an attempt is made to take all these effects into 
account. A dynamic study employing the powerful new tool of the "finite 
element method" (FEM) is clearly the best practical way to perform the 
analysis and such an investigation is currently under way as a thesis study 
by LT. H. L. Crego, USN. Lacking the results of such a study, a "scrambling 
effort" is made in Appendix C to provide a reasonably accurate method of 
estimating natural frequency, accounting for non-uniformity of section and 
for the several non-elementary mechanisms which tend to depress the 
frequency. Briefly, the recommendation in Appendix C is that the frequency 
be calculated for the non-uniform cantilever by use of elementary theory 
(some curves presenting the results of such calculation are included, cf. 
Figure C-5 in Appendix C) and that the depressing mechanisms be accounted 
for (approximately) by use of a frequency reduction factor. 
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This brings us to the third of the subproblems listed above. In his 
original study, Murdock (6) made the requirement r < 0.8, where r = ratio 
of excitation frequency to natural response frequency, basing this recom- 
mendation upon a discussion he had on this subject with Professor J. P. 
den Hartog. Because of the uncertainties surrounding the calculation of 
these frequencies, particularly the effects of rotatory inertia, shear 
deflection, and foundation compliance and inertia, the writer feld strongly 
that the requirement r < 0.8 should be reduced to r < 0.4 and he so recom- 
mended in his original report (3). In the supplement thereto (4) the 
writer made a first attempt at including some of the depressing effects, 
and, as a result, proposed raising the limiting value of r to 0.65. In 
the present study, reported in Appendix C hereof, the completeness of the 
estimate of the depressing effects is believed to be much better. With 
the degree of uncertainty markedly reduced, it is reasonable to liberalize 
the requirement on the ratio r to its original value, namely, again we 
require that r < 0.8. However, it is very important to note that the 
denominator in the expression for r must adequately account for the depress- 
ing effects of rotatory inertia, shear deflection, and foundation compliance 
and inertia. An appropriate way of doing so is by use of the frequency 
reduction factor. Finite element studies, currently in progress, should 
permit refining the analysis. 

The pressure criterion and the non-resonance criterion are the most 
important criteria. However, we are also concerned with the gross effect 
and the fatigue effect of bending. Although the bending moment is obviously 
greatest at the root, for a tapered section, the section modulus varies in 
such a way that the maximum bending stress may not occur at the root as 
was assumed by Murdock (6) and in the writer's earlier analysis (3), the 



10 



latter, incidentally being marred by an analytical error. In Appendix D 
the matter of maximum bending stress and maximum stress intensity is inves- 
tigated. It is found that the previous assumption that the maximum occurs 
at the root is indeed true except for thermowells which are sharply tapered 
or strongly shielded or both. Appendix E similarly studies the fatigue 
effects; in the analysis of fatigue stresses a stress intensification fac- 
tor of 6.0 is assumed. This value is quite conjectural. It is believed 
to be conservative. However the opinion of engineers who deal daily with 
stress intensification factors is definitely solicited on this choice. If 
a value different from 6.0 should be recommended by competent authority, 
the numerical factors in the formula in Appendix E should be proportionately 
modified. 

Appendices D and E represent no changes in philosophy as compared to 
their first presentation in the writer 1 s earlier report. However, the new 
presentations include the effect of partial shielding from the fluid stream, 
something not taken into consideration earlier, and they correct a simple 
algebraic error which was introduced earlier and which was "incorrectly 
corrected" in some subsequent correspondence with the PB51 Committee. 

The bulk of the analysis is presented in Appendices A through G 
hereof. Appendices A through E have been referred to above. Appendix F 
is a Plan and Sequence of Calculations, showing all the recommended 
criteria. Appendix G presents a Numerical Example. 

The writer joins others who may object that the present study, which 
is almost exclusively theoretical, fails to reflect operational experience. 
The most earnest solicitation of information dealing with failure, leakage, 
malperformance, etc., of thermowells attributable to mechanical/structural 
considerations turned up only one case, that of a thermowell which began 
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leaking around at the root connection after many years of successful 
service; in this case, the difficulty seemed almost certainly attributable 
to a poorly executed attaching weld. Mr. J. E. Leary has expressed the 
hope that lack of reports of other cases may indicate a corresponding lack 
of operational difficulties; the writer’s past experience does not lead to 
as sanguine a hope, only to the conclusion that a previously observed 
reluctance to report or to discuss or even to reveal failures of any kind 
of any equipment in any service extends also to thermowells. However, the 
one instance cited above permits the writer to discourse upon his very 
strong conviction that thermowells, as well as any other devices or 
appurtenances the installation of which involves penetrating the pressure 
boundary, must be attached by full penetration welds (or their "equivalent" 
whatever that may be in a particular situation) and that, furthermore, 
the "branch connection reinforcement rules" must be satisfied. See the 
discussion in Section 11 of Appendix A, hereof. 
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1. Basic Analysis 

We consider the elastic-plastic 
behavior of an externally pressurized 
hollow circular cylinder of inner 
radius a and outer radius b. We 
take this to be a one-dimensional 
problem with principal stresses 
cjg, varying only with r. (There is 
no axial or circumferential varia- 
tion.) We assume an ideal elastic- 
plastic material which satisfies 
the usual equations of linear elasticity in the "elastic region" and a 
generalized Tresca condition (Guest's theory) in the "plastic region." 

We assume, in general, external pressurization P sufficient to cause 
plastic behavior for a < r < c and elastic behavior for c < r < b. 

The most fundamental relation which must be satisfied is that of 
radial equilibrium. 

r(da /dr) = a - a (1) 

r 0 r 

which is easily derived by use of a free body bounded by two normal planes 
having unit separation, two cylinders having radius r and r + dr respectively 
and two planes containing the axis and at a dehedral angle of d0 with each 
other. 

We assume that the radial pressure is q = - (a ) at the interface 

r r=c 

between the elastic and plastic regions. In the plastic region we assume 

(this can actually be shown) that cr < a < o <0. Thus the Guest or 

0 z r 

Tresca condition is 



EXTERNAL 




Fig. A. 1 Section of partially plastic 
cylinder under external pressure. 
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o - 0 Q = o 

r 6 



-( 2 ) 



where o is the yield stress for the ideal elastic plastic material. The 

solution which satisfies (1) and (2) and the boundary condition (c ) = 0 

r r — a 

is 



o = -a In (r/a) ; a. = -a[l + ln(r/a)] 
r 0 



(3a, b) 



This gives 

q = Gin (c/a) 

as the interface pressure and it also gives 

(° A ) ^ - -q - a 

0 r— c 



(4) 

(5) 



Lame's solution for the elastic region incorporates the usual elastic 
stress strain relations and the equilibrium relation (1) . It gives 



99 2 9 2 2 2 

= [qc -Pd + (P-q) b c / r ] / (b -c ) 

9 9 9 2 2 2 2 

G 0 = [qc -Pb - (P-q) b c /r j / (b -c ) 

From (6a) it is easily verified that (a ) = -q and (a ) 

r r=c r 

(6b) we find 



r=b 



(6a) 

(6b) 

-P. From 



(0 o ) = [q (b 2 +c 2 ) -2Pb 2 ] / (b 2 -c 2 ) 

0 r=c 



(7) 



and equating this to the expression given by (5) , one gets 

P = q + 0(b 2 -c 2 )/2b 2 (8) 

Using (4) one can also write 

P = o [In (c/a) + (b 2 -c 2 ) /2b 2 ] (9) 

and 

P = q + ( 0 / 2 ) [1 - (a 2 /b 2 ) In (2o/q) ] (10) 

Equations (9) and (10) give P explicitely in terms of c and q 
respectively. However, given P, a more difficult evaluation is required 
to find c and/or q. 
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2. Elastic Pressure P and Ultimate Pressure P 

E U 



If P is sufficiently small, c = a and q = 0. The Lame solutions 

become 



(b 2 -a 2 )a 

r 


= -Pb 2 (l-a 2 /r 2 ) 


(lla) 


(b 2 -a 2 )a 0 


= -Pb 2 (l+a 2 /r 2 ) 


(lib) 


(b 2 -a 2 ) o 


II 

1 

s 

ro 


(11c) 



z 



The third of these is obtained by assuming an end of the cylinder is closed 
and that the external pressure also acts on the closure. Obviously the 
maximum principal stress difference is 

a - a = 2Pa 2 b 2 /r 2 (b 2 -a 2 ) (12) 

r 0 

which is a maximum at r = a, viz. 

(a -a.) = 2Pb 2 / (b 2 -a 2 ) (13) 

r 0 r=a 

Thus, as P is increased, the Tresca condition is first encountered when 

P = P = o(b 2 -a 2 )/2b 2 (14) 

E 

The subscript E indicates elastic action for P < P . 

On the other hand, for sufficiently high pressurization, the inter- 
face radius c assumes the value b and the interface pressure q assumes the 
value P, whence 

P = P y = aln(b/a) (15) 

3. Dimensional Changes 

The subscript U denotes ultimate . However, in contradistinction 
with the usual case to which this word is applied, in the present case 
application of the "ultimate" load does not imply "plastic collapse" with 
"large" deformations. The reason for this is that the geometrical changes 
due to external pressurization are such as to increase the wall thickness. 
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Bridgman (3) has made a large strain analysis of externally pressurized 
cylinders and reference will later be made to his analysis. For the present, 
however, the following simple analysis will suffice. 

Experiments by Bridgman and others and analysis by Bridgman indicates 
vanishingly small axial deformations even under pressures which result in 

gross change of diameters. Under plastic action, the volume remains con- 
stant. Thus, assuming initial radii a , b , and fully plastic action, new 

radii a , b are obtained such that 

2 2 

P = aln(b /a ); b 2 -a 2 = b 2 -a 2 (16a, b) 

2 2 11 2 2 

The first of these reflects fully plastic action (with no strain hardening) 
and the second reflects the volume constancy. We have assumed P > = 

aln(b^/a^), and, given P, a, a^, and b^ we wish to calculate a o and b^ . 

Using the notations 

n = P/Py > 1; a = a 1 / b 1 < 1 (17a, b) 

we easily find 

b 2 = b / (1-a 2 )/ (l-a 2n ) ; a 2 = b^ 11 (18a, b) 

For example, with a^ = .3, b^ = 1.0/ a pressure three times as great 
as P^ gives d = .3, n = 3. Using equations (18) we find b^ = .9543, a^ = 
.0258". The internal radius has been made quite small but equilibrium has 
been restored. If strain hardening occurs, the effective value of n is 
reduced and the distortion is not quite as great as indicated. 

There is nothing, except limitations of pressurization facilities 
to restrict the value of n. For example, assuming a = 65000 psi (accounting 
for heat treatment and some strain hardening) and an applied pressure 
P = 400000 psi, the value of P^ is 78260 psi so that n = 5.11. Then we 
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calculate = .9539", a^ = .0020". Clearly the final dimensions are not 
particularly sensitive to the value of n if n > 2. This calculation is 
consistent with experiments reported by Bridgman (3) except that Bridgman 
observed some cases for which the central cavity closed completely. Obviously 
it would take only a very, very small longitudinal contraction to cause this 
to happen. 

The point of these recent remarks is to the effect that volume con- 
stancy acts to provide dimensional changes which restore equilibrium in 
the case of external pressurization whereas, in the case of internal pres- 
surization, volume constancy gives a thinning of the wall which, in itself, 
acts to remove the situation even farther from equilibrium which can be 
restored, if at all, only by virtue of strain hardening. Thus catastrophic 
collapse jls_ possible for internal pressure but it is not possible for exter- 
nal pressurization. 

Accordingly, simply from the standpoint of maintaining pressure 
integrity, there is no theoretical limit to the external pressure which 
may be applied. However we are also concerned with maintaining a reasonable 
approximation to the original internal dimensions so that the thermocouple 
assembly may be withdrawn and replaced even when the pressure is applied. 

For this reason, we now consider the dimensional changes which can be 
expected with P = P . The exterior surface is barely at the plastic stage 
so that we can use elastic formulas. Experiment indicates that we should 

take e =0, along with a = -P. a n = -a + a = -a (1 + lnb/a) . We find 

z r U 0 r 



0 = Eg 



a 




z 



z 




= - ( 1-v 2 ) a + (l-v-2v 2 ) a 



r 



= - (1-v 2 ) a - ( 1-v- 2 v 2 ) a lnb/a 



(19) 
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These values obtain at r = b with P = P^. Thus the radial deformation at 
r = b is 

6b = be = - [l-v 2 + (l-v-2v 2 ) Inb/a] (20) 

0 E 

There is an ambiguity concerning the value of v to employ. For P slightly 
less than P^, v=. 3 (approximately) whereas with P slightly greater than P y , 
v=.5 but the elastic equations are not applicable. However, for our present 
purposes, the effect of v does not alter the conclusion at which we shall 
arrive. We have 

r -u. 

e 0 = — = -(a/E) ( . 91+ . 521nb/a) (v=.3) 

= - . 75a/E (v=. 5) 

- 3 

Taking, for example, a reasonable value of a/E = 1.2 x 10 and b/a = 3 
we have 

— = -.0018 (v= . 3) 

b 

= -.0009 (v= . 5) 

The larger value is probably more accurate and for our present purposes it 
is conservative. 

Assuming volume consistency we have 

2 2 2 2 

(b+6b) - (a+6a) = b - a (21) 

so that, to first order 

6 a/a = (b/a 2 ) 6b = (b/a) 2 6b/b = -.0162 

Thus, in a typical case, if a = .16", and b = .48", then 6a = -.026". That 
is the interior radius shrinks by .026". This order of magnitude appears 
tolerable and our conclusion is that if the applied external pressure does 
not exceed P^ then the decrease in internal diameter will not adversely 
affect the ability to remove and replace thermocouple elements. 



20 



4. Repeated Pressurization, Elastic Shakedown, and the One-Cycle Shake- 
down Pressure P*. 

We assume initial pressurization to a pressure P (P < P < P ) so 

E U 

that there has been plastic behavior from r = a to r = c, (a < c < b) . We 
presume however that subsequent removal of this pressure results in no 
additional plastic behavior. That is, the depressurization operation is 
purely elastic. Thus we can arrive at the final state of stress by super- 
posing the Lame stress system 



= Pb 2 (l-a 2 /r 2 ) 


l / (b 2 -a 2 ) 


(22a) 


= Pb 2 (l+a 2 /r 2 ) 


1 / (b 2 -a 2 ) 


(22b) 



upon the system given by equations 3 and 6. The maximum (equivalent or 
Tresca) stress state after depressurization occurs at r = a and is positive, 
i.e., at r = a we have 



a = 0, a = -a + 2Pb 2 /(b 2 -a 2 ) 
r o 



(23) 



If the yield condition is not to be exceeded (this time in the 
opposite sense from originally) , we must have 



(a ) - (a ) < a 

0 r=a r r=a 



(24) 



and this gives 



Pb 2 /(b 2 -a 2 ) < 5 ; P < P* = 2P = a (b 2 -a 2 )/b 2 

E 



(25) 



We will use the symbol P* = 2P^ and refer to the condition described 
above as "one cycle shakedown" since it assures that after the plastic 
yielding occurring on initial pressurization there can be no subsequent 
yielding. 

Our previous analysis has led us to conclude that if P < P^ the 
deformations will be acceptable. Thus, we obviously wish to compare the 
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pressures and P*. Equating these values, and using the notation a = a/b, 
we immediately obtain the equation 

1 + lna = a 2 (26) 

which has two roots, a = 1, which is meaningless for our application, and 

a = 0.4503 (27) 

For a > 0.4503, P^ < P* and if we required P < P^, then one cycle shake- 
down is absolutely assured. For a < 0.4503, one cycle shakedown requires 
limiting P to a maximum of P* < P^. 

5. Basis of Recommended Criteria 

We can represent our principal findings in a very compact form. Using 
the notation a = a/b, we have 

P*/o = (1-a 2 ) , P^/o = -ln(a) (28 a,b) 

Our recommendation , 

P £ Min{P^ J , P*} = Min{-ln(a) , (1-a 2 )} (29) 

is indicated by the heavy line in Fig. A-2. This condition assures both 

(a) not exceeding the ultimate pressure P , so that we are assured 
of tolerable dimension changes under pressurization, and 

(b) one cycle shakedown so that all yielding occurs on initial 
pressurization and none on subsequent depressurization and 
repressurization. 

This recommendation appears reasonable at this time. Since a < .4503 
for current designs of thermowells for high pressure applications, the 
criterion is conservative since surely there is no compelling reason to 
call for one cycle shakedown. Safe operation would also be assured if 
shakedown were to occur in two, three, or any other relatively small number 
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2.0 




Radius ratio a = a/b 

Fig. A- 2 Summary of most significant pressure calculations. 

of cycles. P/a curves for such n-cycle shakedown have been con jecturally 
sketched in Fig. A-2 for n = 2,3, and 4; these curves should not be used 
quantitatively. The analytical difficulties involved do not presently 
warrant working out their actual shape. However, in the future, in any 
case where the presently recommended criterion should prove to be restric- 
tive there would be ample reason to reconsider the matter so as to permit 
two or three cycle shakedown or so as to permit exceeding the pressure P . 
6. Numerical example 

So as to demonstrate the self consistency of the preceding analysis, 
it is desired to consider a practical case. We take a/b = a = 0.3 and 
easily calculate P E = 0.4550a, P* = 0.9100a, and P^ = 1.2040a. We take 
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P = p* and let 3 = c/b. From equation 9, which becomes 



l-a 2 = In ( 0/<x) + (l-g 2 )/2 (30) 

we calculate 3 = .51655. This also gives q/a = InB/a = .5434. Stress 
calculations are shown in Table A-l and Figure A- 3. The "final state" 
referred to is at the end of the first pressurization cycle, external 
pressure P = .9/a having been applied once and then removed. After initial 
yielding, subsequent application and removal of P = .9/0 ceases stresses 
to vary between solid and dashed extremes. 



TABLE A-l 



Stress Calculations for Numerical Example 





PRESSURIZATION 


REMOVAL 




FINAL 


STATE 


Radius 


V” 


o /o 
r 


Q 

CD 

1 

Q 
Q 1 


V 5 


Q 

\ 

Q 1 


V 5 


a /a 
r 


(o.-a )/o 
0 r 


.3b 


-1.000 


0 


-1.000 


2.000 


0 


1.000 


0 


1. 000 


.4b 


-1.288 


-.288 


-1.000 


1.562 


.438 


.274 


.150 


.124 


.5b 


-1.511 


-.511 


-1.000 


1.360 


.640 


-.151 


.129 


-.280 


c 


-1.543 


-.543 


-1.000 


1. 337 


.663 


-.206 


. 120 


-.326 


.6b 


-1.414 


-.673 


- .741 


1.250 


.750 


-.164 


.077 


-.241 


. 7b 


-1.316 


-.771 


- .545 


1.184 


.816 


-.132 


.045 


-.177 


.8b 


-1.252 


-.835 


- .417 


1.141 


.859 


-.111 


.024 


-.135 


.9b 


-1.208 


-.879 


- .329 


1.111 


.889 


-.097 


.010 


-.107 


b 


-1.177 


-.910 


- .267 


1.090 


.910 


-.087 


0 


-.087 



7. Remarks About Experiments and Other Analyses. 

In a previous analysis (4) , the writer arrived at identical results 
and conclusions but the details may appear to be different than those given 
above. The reason is that in a misguided effort to make the presentation 
in terms of familiar material, the previous analysis (4) dealt with interior 
rather than external pressurization. Since the "failure law" (i.e., the 
Tresca condition) is independent of the hydrostatic stress state, i.e., is 
independent of the first scalar invariant of the stress tensor, internal 

pressurization is equivalent to external depressurization, and, with an 
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Fig. A- 3 Stress distributions in numerical example. 

overall change in algebraic sign, the latter is equivalent to external 
pressurization . The fact that the same conclusions are reached on the 
basis of the present analysis as were reached in (4) is sufficient evidence 
of the equivalence. However, instead of finding this viewpoint simplifying 
or illuminating, many readers of the earlier report were dismayed and con- 
fused. Accordingly, in the present analysis the actual, rather than an 
equivalent situation is treated. 
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The procedure here employs the simplest and most common plastic analy- 



sis to be found in the literature. The use of the Guest or Tresca condition 
is consistent with usage in the A.S.M.E. Boiler and Pressure Vessel Code. 

A readily available development (of the interior pressurization problem) 
is to be found in Timoshenko's popular textbook (12). The analysis is 
probably originally due to Nadai (9). However, it is not the only view- 
point. Bridgman (3) develops a generalized version and incorporates large- 
strain analysis; his development is essentially for fully plastic behavior. 
Hill et al. (6) consider modifications of the present analysis based on 
resolution of an undesirable discontinuity of axial strain at the elastic- 
plastic interface; their analysis predicts slightly less radial deformation 
than does the analysis given here for a given pressurization. A later work 
by Nadai (10) devotes three chapters to analysis of pressurized cylinders. 

A very recent work by Save and Massonnet (11) summarizes work to date and 
offers a large bibliography; they cite additional recent studies. The 
analysis summarized in (11) is the same as that given here. 

Burst tests, such as described by Faupel (5) simply do not apply to 
the problem in which we are interested. However such tests have served to 
verify the general reliability of all the plastic analyses available; strain 
hardening is such as to mask differences. 

Hill et al. (7) provide an analysis which indicates that if the 
Mises rather than the Tresca condition is used, is increased by 15%, 
i.e., P^ = (2//3) aln (b/a) . Thus, use of the Tresca condition appears to 
be conservative. 

Throughout the analysis and discussion to this point we have assumed 
axial symmetry. Specifically, we have not considered a mode of collapse 
in which the section becomes ovalized or goes out-of-round. The ASME 
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"rules" for externally pressurized (thin wall) circular cylinders are based 
upon the predication of out-of- round deformation. Timoshenko and Gere (14) 
describe the genesis of the ASME procedure as combining a classical shell 
buckling analysis with what essentially amounts to the present formula (15) , 
which, for thin wall, takes the form 

P rT = a In [ (l+t/2r )/(l-t/2r )] t ta/r (31) 

U m m m 

as is given by the most elementary analysis. Ref (14) indicates also that 
the ASME rules of 1933 used an artif ideally low value of a = 26000 psi; 
presumably modern versions of the ASME rules, which now provide curves 
for a number of different metals of engineering importance, reflect more 
realistic values of a. However, there do not appear to be any analyses in 
the literature which treat out-of-round buckling of externally pressurized 
cylinders having ratios a = a/b as large as those employed in thermowells 
for high pressure service. 

Confining attention to thermowells for high pressure service, there 
seems to be no "engineering sense" in applying any criterion more restric- 
tive than the one recommended in the preceding paragraphs. Externally 
pressurized cylinders having thick walls simply do not go out-of- round. 
Bridgman (3) conducted a number of tests on tubes of various steels having 
O.D. = .3125, I.D. = .0998. This corresponds to a = .32 which is approxi- 
mately the value used in the examples in the present analysis and also 
corresponds to current industrial practice for high pressure installations. 
He subjected these tubes to as high as 412000 psi external pressure. In 
each case the tube simply decreased in diameter, while maintaining its 
length almost exactly without change. Equation (16b) was satisfied; that 
is, there was no volume change that could be detected. In one case, of a 
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very soft steel under 412000 psi, the central cavity appeared to close up 
completely. However, there was no failure or loss of pressure integrity. 
Thus, for such thick wall cylinders, non symmetric distortion simply does 
not take place. The criterion recommended here (Equation 29) should assure 
that dimensional changes remain acceptably small and that shakedown to 
elastic conditions occurs promptly so that there is no danger of ratchetting 
or low-cycle fatigue. 

8. Thin Wall Thermowells 

The basic motivation behind this re-examination of design criteria 
for thermowells resides in the integrity of thermowells against very high 
external pressures, and the analysis thus far has been of such cases. 
However, design rules should cover all possible pressurizations, including, 
as an example, thermowells in exhaust gas ducts. However, it is not within 
the scope of the present study to consider such cases seriously or in 
detail. It seems reasonable to suppose that for case with "nominal" 
values of external pressure, the current interpretation of the Power Test 
Code, whatever that may be, should be considered applicable. In the appro- 
priate part of this document (1) the maximum gage pressure is given by a 
formula P = K^S where Kj is a constant varying between 0.155 (for "large" 
thermocouple elements) to 0.412 (for "small" elements). If this criterion 
is satisfied for a thermowell having (roughly, say) the dimensions given 
in this document (1) , there seems to be no reason to question the accepta- 
bility of the design on the basis of pressurization. For designs which 
must vary from the dimensions indicated in (1) , it seems reasonable to 
attempt to apply the rules in the Unfired Pressure Vessel Code (2) . 

The only questions which appear to remain (on the question of wall 
thickness vs. pressure) are these. (1) For very low external pressuriza- 
tion, a certain degree of structural strength, possibly more than might 
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be required by other criteria for thermowells, might be required to with- 
stand the loads applied during shipment and installation or due to inadver- 
tently applied mechanical loads during operation and maintenance, and (2) 
for "medium" pressure situations under what circumstances should one be 
concerned about true buckling in which ovalization or lobar deformation 
occurs . 

We shall not concern ourselves further with (1) above. Let us turn 

attention to (2) . The smallest value of D /t contemplated in the ASME 

o 

Unfired Pressure Vessel Code (2) (Appendix V) is D /t = 10, and as was 

o 

pointed out in (14) , the criterion in this case is gross yielding with no 

change from the circular shape, using a substantial (>2) factor of safety. 

In the writer* s opinion this is greatly overconservative for thermowells. 

However, let us not argue with it. We propose therefore that for D /t ^ 

o 

2 

10, the Unfired Pressure Vessels be employed. The parameter D /t = 

o 1-a 

using the parameter a = a/b introduced earlier. Thus D /t = 10 corresponds 

o 

to a = .8. We believe that the strict application of the UFPV rules for 

D /t < 10 may be uneconomically conservative. Accordingly, we suggest that 
o 

the entire gamut of pressure criteria be based as follows: (1) rules pre- 
viously suggested for 0 = a = .6; (2) UFPV Code rules for .8 = a = 1.0, 

(3) linear interpolation between. 

9. Factor of Safety 

Although the preceding discussion indicates no need for a "factor of 
safety" since catastrophic dimension change is not possible and even con- 
siderable overpressurization can result in nothing worse than squeezing 
down on the thermocouple element, nevertheless it is customary to provide 
for a factor of safety or its equivalant to account for inadvertent occa- 
sional overheating and/or overpressurization, the possibility of individual 
metallic specimens failing to possess the physical properties called for in 
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